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Abstract 

We compare two existing approaches to calculating the decay of molecular quarkonium states 
to conventional quarkonia in effective field theory, using X(3872) —)• XcJT^^ as an example. In one 
approach the decay of the molecular quarkonium proceeds through a triangle diagram with charmed 
mesons in the loop. We argue this approach predicts excessively large rates for r[X(3872) —)■ XcJT^^] 
unless both charged and neutral mesons are included and a cancellation between these contributions 
is arranged to suppress the decay rates. This cancellation occurs naturally if the V(3872) is 
primarily in the 1 = 0 DD* + c.c. scattering channel. The factorization approach to molecular 
decays calculates the rates in terms of tree-level transitions for the D mesons in the X(3872) to 
the final state, multiplied by unknown matrix elements. We show that this approach is equivalent 
to hadronic loops approach if the cutoff on the loop integrations is taken to be a few hundred MeV 
or smaller, as is appropriate when the charged D mesons have been integrated out of the effective 
theory. 
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I. INTRODUCTION 


The last ten years have seen a plethora of discoveries of unconventional quarkonia, ^ 
the hrst and most studied of these being the X(3872) [I]-|l]. Because of its proximity 
to the D^D*^ threshold it is thought by many authors to be a molecular state. If the 
state consists primarily of the C even linear combination of neutral D mesons, D^D*^ 
+ C.C., the binding energy is —0.11 ± 0.21 MeV, and this state is a very shallow bound 
state. For the central value of this binding energy, one calculates the typical separation of 
the and D*^ to be approximately 10 fm, which is an astonishingly large length scale 
compared to typical hadronic scales. Ref. [S] exploited this separation of scales to construct 
an effective held theory for the X(3872) called XEFT. Heavy hadron chiral perturbation 
theory (HH^PT) [T0l - [T2] is matched onto a non-relativistic theory of neutral D mesons and 
pious. Their interactions are constrained by the heavy quark and chiral symmetries of QCD. 
A contact interaction is tuned to produce a shallow bound state in the D^D*^ + c.c. channel 
which is the X(3872). The structure of the theory is similar to effective held theories of the 
deuteron and low energy two-body nuclear physics HSIE]. 

For processes that are dominated by long-distance aspects of the X(3872), such as 
X(3872) —or D^D^'y, this theory reproduces ehective range theory (ERT) at low¬ 
est order. ERT predictions for these X(3872) decays were hrst calculated in Refs. [T^ IT6] . 
XEFT allows for the systematic inclusion of corrections to these predictions from pion loops 
and higher dimension operators. Ref. [S] showed the corrections from pion loops were neg¬ 
ligible, at least for the process X(3872) —)■ D^D^ir^. The ehect of hnal state interactions on 
the reaction X(3872) —was recently studied in Ref. [T7j. For calculations of many 
processes within XEFT, see Refs. XEFT has also been used to calculate the quark 

mass dependence of the X(3872) binding energy in Ref. [25], for a related EFT calculation 
see Ref. |2S]- 

Many observations of X(3872) involve decays to conventional charmonia, including 
X(3872) —>■ J/'ip7T~^7r~, J/'tpTT~^7i~TT^, '7/V’7, and 'ip{2S)'y. The X(3872) has also recently been 
observed in the decay of the exotic quarkonium state F(4260) —)■ X(3872)7 [27] . Ref. [2H| 
predicted an enhanced rate for the decay y(4260) —)• X(3872)7 based on the assumption 

^ For a review of recent developments in quarkonium spectroscopy, we refer the reader to Refs. [SHH] 
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that the Y (4260) is a DDi molecule, while other authors interpret the Y (4260) as a char- 
monium hybrid [22], so this transition probably does not involve a compact cc state. Other 
decay and production processes with conventional charmonia such as J/ip or 'ip{2S) involve 
short-distance scales since the D and D* must coalesce to couple to a conventional char- 
monium. For these decays there exist two distinct approaches to applying XEFT in the 
literature. The approach Erst taken in Ref. [IB] is to use HHyPT to calculate the transition 
of D^D*^ to the hnal state, then match the resulting amplitudes onto XEFT operators. The 
resulting prediction for the partial decay width of the X(3872) is given by an expression of 
the form 


F[X(3872) ^ F.S] oc Oxeft x a[D^D*^ + c.c. ^ F.S .], (1) 

where F.S. denotes the hnal state (which includes a charmonium) and Oxeft is an XEFT 
operator. This operator plays the same role as the wave function at the origin squared 
in a traditional approach to bound state calculations. The numerical value of the XEFT 
operator is unknown and must be extracted from data. Since the and D*^ must coalesce 
to form the compact charmonium, part of the process involves short-distance physics that 
is not determined by the universal nature of the long-distance part of the X(3872) wave 
function, and this physics is encoded in Oxeft- Similar factorization theorems for X(3872) 
decay and production were developed in Refs. janiEi]. We will refer to the approach to 
X(3872) decays advocated in Refs. [IHl ED] El] which yields a factorized formulae of the 
form of Eq. Q as the factorization approach to X(3872) decays. 

The second EFT approach to X(3872) production and decays is advocated in, e.g.. 
Refs. [2H1 E2]- The decay involving the conventional quarkonium proceeds through a 
loop diagram in which both the X(3872) and the conventional quarkonium couple to 
heavy mesons. In this case the X(3872) coupling to heavy mesons in the loop is hxed 
by the residue of the pole in the T-matrix. In some cases [28] a power counting argu¬ 
ment shows that the hadronic loop is lower order than any tree-level diagram and the 
hadronic loop approach is more predictive than factorization since there is no undeter¬ 
mined XEFT matrix element. Whether or not this happens depends on the quantum 
numbers of the states involved in the transition. For example, in the radiative transitions 
J/'ip,'ip{2S) —)■ X(3872)7 a counterterm appears at leading order, so it is not possible to 
predict the ratio F['0(2S') —)■ X(3872)7]/F[J/'0 —)■ X(3872)7] [32]. A similar conclusion was 
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reached in the factorization approach in Ref. [21] . Even when there is no tree-level countert¬ 
erm at leading order, this approach still may not be entirely predictive since the couplings 
of heavy mesons to conventional quarkonia in the loop could be unknown. We will refer to 
the approach to X(3872) decays in which the decay is assumed to go through a hadronic 
loop as the hadronic loop approach to calculating X(3872) decays. 

In addition to two different approaches to calculating X(3872) decays, there are also dif¬ 
ferent choices of the relevant degrees of freedom appropriate for an effective theory suitable 
for describing the X(3872). In the literature there are calculations within both the factoriza¬ 
tion approach and the hadronic loop approach that only include neutral D mesons as explicit 
degrees of freedom, since the X(3872) is considered a shallow bound state of these mesons 
alone. ^ Refs. [^IM] have emphasized the importance of including charged D mesons as well 
in the calculations of the decays X(3872) —)■ J/V'T, , and J/'07r+7r“7r°. Note that 

the charged meson threshold D~^D*~ + c.c. is considerably farther away from the X(3872) 
mass than the neutral threshold. The binding energy is 8.2 MeV and the corresponding 
estimate of the separation of the charged mesons in the X(3872) is 1.1 fm. This is roughly 
a factor of 10 smaller than the central value for the corresponding estimate for the neutral 
channel. The charged D mesons are separated by a distance that is not much larger than the 
size of the hadrons themselves. At length scales larger than 1 fm, the wavefunction is cer¬ 
tainly dominated by the neutral mesons. In the original formulation of XEFT the charged 
mesons are integrated out of the theory and their effects subsumed in to short-distance 
XEFT operators. However, for processes like decays to conventional charmonium, in which 
both long and short distance scales are important, it may be desirable to include these as 
explicit degrees of freedom. 

The purpose of this paper is to compare the different approaches to calculating the 
decay of X(3872) to conventional quarkonia, using the decays X(3872) —>■ Xcjtt® as an 
example. These decays were hrst studied in Ref. [33] where it was pointed out that the 
relative rates for different J are predicted by heavy quark symmetry and can be used to 
distinguish between different interpretations of the X(3872). These decays were studied 

^ For some processes there are other justifications for neglecting the charged mesons. For example, in 
the calculation of X{3872) radiative decays in Ref. [25], the charged mesons were neglected because the 
neutral charmed mesons couple much more strongly to the photon. However, charged mesons are included 
in analysis of 2f(3872) radiative decays in Ref. [52] . 
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in XEFT in Ref. [18], which used factorization in a theory with only neutral D mesons as 
explicit degrees of freedom. Ref. [TH] showed that within this approach there are two distinct 
long-distance and short-distance mechanisms contributing to the decay and the relative rates 
depend on the relative importance of the two mechanisms. The authors of Ref. [IH| also 
computed the partial widths using the hadronic loop formalism, with only neutral D mesons 
as explicit degrees of freedom, but as we will see in the next section, this yields exceedingly 
large partial widths for X(3872) —)■ that are in conflict with experiment, so this 

approach was discarded and the calculation was not published in Ref. [IH]. This result 
is somewhat model dependent as the predicted rates depend on the unknown coupling of 
the XcJ to charmed mesons, which is estimated using the model in Ref. [36] • However, to 
make the predicted partial widths for X(3872) —)■ consistent with experiment requires 

that this coupling be almost two orders of magnitude smaller than what one expects from 
naive dimensional analysis. We conclude that the hadronic loop approach with only neutral 
charmed mesons as explicit degrees of freedom is inconsistent with experiment. The hadronic 
loop approach can be made consistent with data if charged mesons are included as explicit 
degrees of freedom. If the X(3872) has nearly equal couplings to the charged and neutral 
channels a cancellation between charged and neutral loop diagrams suppresses the rate. This 
cancellation occurs naturally if the X(3872) is an / = 0 state. An J = 0 interpretation of 
the X(3872) has been put forth by other authors [33l|3l| and is consistent with the observed 
ratio r[X(3872) —)■ J/'^7r+7r“7r°]/r[X(3872) —)■ = 0.8 ± 0.3 [2Z| if one accounts 

for differences in two- and three-body phase space [35] . 

In section III, we discuss how the hadronic loops approach is related to the factorization 
approach. We show that the hadronic loop integral can be expressed as the convolution of the 
ERT wave function of the X(3872) with the tree-level matrix element for D*^D^ —)■ Xcjtt® 
We simplify the calculation by dropping some terms 0{p^/{mr)ET,)) which only changes 
answers by a few percent. Then the hadronic loop integrals contain contributions from 
two very different scales, 7 ^ = 14 MeV and \/rrii)E^ 850 MeV. Here jn is the binding 
momentum in the neutral channel, niD is the D meson mass, and is the energy of the pion 
in the decay. The contribution from loop momentum of order gives the dominant 

contribution to the integral, but we argue that in a theory in which the charged mesons 

^ Here and throughout this paper charge conjugate channels are implied. 
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have been integrated out, the theory must be thought of as having a cutoff A ~ 7 c, where 
7 c is the binding momentum in the charged channel, since the ERT form of the X(3872) 
wave function, with only neutral D mesons, is no longer reliable above this momenta. If 
the hadronic loop integral is performed with a cutoff. A, such that 7 „ <C A <C y/moK^, 
one recovers the results from the factorization formalism. We also show that if the theory 
contains both charged and neutral D mesons and the cutoff is taken to be large compared to 
and the couplings of the X(3872) to the charged and neutral channels are equal, 
then the large contributions from the part of the hadronic loop integrals cancel 

and the remainder is well approximated by the factorization formulae, with A ^ 7r7c/2 200 

MeV. In the hnal section we give our conclusions. 

Our study is closely related to that in Ref. [39] which compared the wavefunction at the 
origin squared prescription to the hadronic loop approach in hadronic molecule decays to 
two photons. Their main conclusion, relevant to this paper, is that when the range of the 
forces binding the hadronic molecule is much smaller than the distance scale associated with 
the annihilation, the hadronic loop approach is appropriate, while the wavefunction at the 
origin prescription is appropriate in the opposite limit. This is consistent with our analysis, 
but it is unclear whether the assumptions appropriate to the hadronic loop approach apply 
in the case of the X(3872). The momentum scale characterizing the annihilation process is 
■y/mjjEjr ~ 850 MeV, corresponding to a length scale of ~ 0.23 fm, which is comparable to 
the size of the charmed mesons themselves. It is not clear a priori that the ERT wavefunction 
of X(3872) will be correct down to such a short distance, but if it is then charged charmed 
mesons must be included as explicit degrees of freedom. If the ERT wavefunctions are only 
valid for much longer distance scales, than a factorization approach may be more appropriate. 
Hopefully, future experimental and theoretical studies will clarify which approach is more 
suitable for X(3872). 
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II. HADRONIC LOOPS 


In this section we will consider the X(3872) decays to Xcj in the hadronic loops approach. 
The LO HHyPT lagrangian for the charmed mesons is 

+ Ai[Hl{iDo)abH,] + g Ai[HI a ■ + ^Tr[iyt a* H, a']. (2) 

We use the two component notation of Ref. [ID] . The held Ha is given by 

Ha = Va ■ a + Pa , (3) 

where I 4 annihilates D* mesons and Pa annihilates Da mesons. The subscript a is an SU (2) 
index, and a = 1 for neutral D mesons. The corresponding held for antimesons is Ha- The 
held Aab is the axial current of chiral perturbation theory, Aab = —VTiab/fw + where 
is the pion decay constant and Tiab are the Goldstone boson helds. The lagrangian coupling 
the Xcj to heavy mesons is 

= i^Tiix^^Haa^ Ha] + ^Tiix^^Haa^ HbhjkA’^ab + h-c., (4) 

where the Xcj helds are represented by 

+ (5) 

The transformation rules for the various helds under the symmetries of the theory can be 
found in Ref. [18] . 

In the hrst part of this section we will include only the neutral D mesons as physical 
degrees of freedom. The hadronic loop diagrams for the decays of X(3872) to the XcJ are 
shown in Fig. In this hgure the black line represents the interpolating held for the X (3872) 
and the gray lines are the Xcj- The internal lines are the neutral D mesons, with a single line 
representing the or and a double line for the D*^ or D*^. For power counting we use 
the V counting of Refs. [2HI |32], which is appropriate for the hadronic loop approach. The 
couplings of the X(3872) and the Xcj to the D mesons have no derivatives, so these scale 
as The pion is derivatively coupled so that interaction scales as In the loops, the 
integration measure scales as and each propagator scales as so the diagrams scale as 
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FIG. 1: Hadronic loop diagrams contributing to the decays X(3872) —>■ XcJT^^ 


Ptt/v- There is also a loop diagram that contains a bubble with the four-particle interaction 
multiplied by ci in Eq. Q. This diagram contains one fewer propagator, the four-particle 
interaction still has a derivative acting on the pion field, so the diagram scales as p-j^v and 
is suppressed by in the v expansion. Finally, there is a possible tree-level X{3S72)-Xcj-'^ 
coupling, which would scale as p-„ and is suppressed by v in the v expansion. Hence, the 
diagrams shown in Fig. [^are the leading contribution to X(3872) —)■ Xcj^^ in the hadronic 
loop approach. 

The coupling of the X(3872) to the D^D*^ + c.c. is 27r"Y^ / p nn* [IT] . where pdd* is the 
reduced mass of the and D*^ and 7 „ is the binding momentum in the neutral channel, i.e., 
7 n = ^/2fJ,J:)D*BEn, where BEn is the binding energy in the neutral channel, m^.o + uido — 
mx( 3872 )- If one uses the interpolating held {D^D*^ + D°D*^)/y/2 to represent the X(3872) 
this factor arises from wave function renormalization obtained using the LSZ formalism for 
composite operators, see, e.g.. Refs. PE2]. Computation of the rates is straightforward and 
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we simply quote the prediction for the rates: 


r[X(3872) ^ 

971 ^;^ mx(3872) 

r[X(3872) ^ X,i7r°] = pi 

DTT^/^ mx(3872) 

r|Jf (3872) ^ X= 2 >r"l = pi F,[',,„ A„, E,]\ (6) 

1877"'/^ mx(3872) 

Here 5 ^ = 0.54 is the axial coupling of the D mesons to the pion, gi is the coupling of the 
Xcj to D mesons, = 130 MeV is the pion decay constant, rn^^j{mx( 3872 )) is the mass of 
the Xcj (-^(3872)), and E^{p^) is the energy (momentum) of the pion in the decay. The 
factors Fi[xn, Etj] come from the loop integration and are given by 

F„|7„,A„,£J = ^f(jI,2^od-{E.-^o) + j,j] (7) 

( 2 2 \ 

( 2 2 \ 

7^,2p,,^*(i7^ + Ao) + ^,^j . 


To simplify Eq. ([^ in some places we have approximated itido Ri m^.o 2hdd*, which is 
accurate to 4%. In Eq. ([^, Aq = rrin-o — rriDO is the hyperhne splitting for the neutral D 
mesons, and the function E{a, b, c) is given by 

1 


E{a,b,c) = 


dx 


\/a + b 


1 

1 

7 c L' 


tan 


-1 


X — cx^ 

b 


sm 


-1 


ac 

b 


tan 


-1 


b — 2c 


\/b‘^ + 4ac 


2y/c\/a + b — c) _ 
b — 2c 


— sm 


-1 


y/b"^ + 4ac 


( 8 ) 


The hrst analytic expression for the evaluation of the integral is appropriate for a,5, c > 0, 
a,c,a + b — c ^ 0, which is always the case for us. In our case we always have b S> a, c, 
and the second analytic expression in Eq. (|^ is better suited for expanding in a/b and/or 
c/b. In the heavy quark limit where the Xcj are degenerate and Aq = 0, Fo[ 7 „, Aq, = 


^ This value for g is obtained using the recent measurement of r[Z?*+] = 83.4 ± 1.8 keVfTTl l43] times the 
measured strong decay branching fractions for the D *+ m\ and the tree-level HHyPT expression for the 
strong decay width of the £)*+. 
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Fi[7n, Ao,-E tt] = -F 2 [ 7 n, Ao,-Ejr], E-^ and are the same for all three decays, and the rates 
are in the ratio ro:ri:r2::4:3:5, where Fj = r[X(3872) —)■ In reality, the 

small hyperhne splittings signihcantly affect the valne of mnltiplying each decay, and the 
factors of p^(m^^_^/mx( 3872 ))-Fj[ 7 n, Aq, 77,^] differ signihcantly between the three decays, so 
we hnd 

Fo : Fi : Fa :: 3.2 : 1.2 : 1.0. (9) 

Ref. [35] calcnlates these ratios by weighting the heavy qnark spin symmetry prediction with 
the p^ factors mnltiplying each decay, obtaining Fq : Fi : Fa :: 4p^ : 3p^ : 5p^ :: 2.7 : 0.95 : 
1.0. The factors of p^ acconnt for most of the deviation from heavy qnark spin symmetry 
predictions, remaining factors give corrections of order 20 — 25%. 

To compnte the absolnte rates in this approach, one needs to know the conpling constant 
Pi in Eq. (|^ and the binding momentnm, 7 „. From the binding energy BEn = 0.11 ± 21 
MeV, we hnd Xn = g MeV. Becanse Xcj is a conventional qnarkoninm rather than 

a bonnd state of charmed mesons the conpling pi is an nnknown parameter. We will nse 
the resnlts of Ref. [36], which estimates the conpling ® by nsing a vector meson dominance 
argnment to hnd gl ~ 'i^xc.ol (^/xco)’ where = (0|cc|xco) and is calcnlated to be 510 MeV 
from QCD snm rnles. Using this estimate, gl = 1/(457MeV), and we hnd 

F[X(3872) ^ Xeovr°] = 3.8 MeV 
F[X(3872) ^ = 1.4 MeV 

F[X(3872) ^ Xe 2 vr°] = 1.2MeV. (10) 

All of these partial widths separately exceed the cnrrent experimental bonnd on the total 
width, Fx < 1.2 MeV [37]. 

The partial widths, F[X(3872) —)■ Xcjti'^], which are presently nnmeasnred, mnst in fact be 
orders of magnitnde smaller than the existing bonnd on the total width. We will next hnd an 
npper bonnd on the snm of the partial widths, F[X(3872) —)■ Xcjt^^]- Theoretical calcnla- 
tions P Hg ini |22] of F[X(3872) ^ T)°Zl°7r°] hnd F[X(3872) ^ = T[D*^ DV] 

in the limit of zero binding energy. Zl°7r°] has not been directly measnred, bnt can 

be obtained using the total width F[Z1*+] = 83.4 ± 1.8 keV [ST] |13| and —)■ Zl+7r°] = 

® Our definition of the coupling gi is a factor of l/-\/2 smaller than the gi defined in Ref. [55] . 
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30.7±0.5% which together give r[Zi)*+ —= 25.6 ±0.69 keV. In the isospin symmetry 
limit, —)■ Zi)+ 7 r°] = —)■ Noting that each decay scales like and taking 

into account differences in phase space, we find —)■ = 36.4±0.98 keV. Therefore, 

we expect r[X(3872) ^ = 36;?o keV. The central value here is our extracted value 

of —)■ which has only a few percent uncertainty from experimental uncertainties 

and isospin violation. The uncertainty in r[X(3872) —is obtained by assuming 
that the binding energy of the X(3872) is between 0 and 0.3 MeV, and using the theoretical 
calculation of r[X(3872) —in Ref. [9], which includes corrections from range cor¬ 
rections, pion loops, and higher dimension operators. Furthermore, the branching fraction 
r[X(3872) ^ D^D%^]/T[X{3872)] > 32% [37], implying r[X(3872)] < 131 keV. (We use 
the largest value of r[X(3872) — D^D^n^] in our quoted range to obtain this bound.) The 
branching ratio for any of the strong decays to conventional quarkonia is considerably smaller 


than this. For example, F[X(3872) — D^D^'k^]/T[X{3872) —] = 8 . 8^35 
plying F[X(3872) ^ = d.ll^-s keV if F[X(3872) ^ D^D^-k^] = 36 keV. The 

total partial width to final states and 'ijj{2S)'j constitute 

at least 39.5% of the total width 1371 , so the total partial width to all other states is less 


im- 


than 79 keV, so X]jr[W(3872) —)■ < 79 keV and we can see from Eq.(lO) that the 

hadronic loops prediction is almost two orders of magnitude too large. If this situation is to 
be fixed by using smaller values of gf and changing no other parameters, we must require 
Qi ± 1/(37 GeV) which seems implausibly small from the point of view of naive dimensional 
analysis. 

One way to fix this is to include both the charged and neutral mesons in the theory, since 
the decay rate is naturally suppressed if the X(3872) couples to charm-anticharm mesons in 
the 1 = 0 channel. Refs. [33l l34] have emphasized the necessity of including both charged 
and neutral D mesons in the context of 77(3872) —)• J/V’ 7 , and J decays. 

When the charged channel is included as well the formulae of Eq. (§ generalize to 

pI (^oT^obn, Ao, E^] - ^+Fo[7c, A+, E^])^(11) 


F[X(3872) ^ x,ovr°] = 




^X(3872) 


F[X(3872) ^ = 

F[X(3872) ^ x, 2 vr°] = 


01 _ 9 ^ 9 l ,4 ^xci 3 


127r3f2 

5g^gf 

3671^ f- 


Pdd* 


^X(3872) 


Pn {9oFi[yn, Ao, E^] - ^+Fi[7c, A+, E„]y 


m 


;PdD* 


Xc2 


^X(3872) 


■pI (^o7^2[7n, Ao, - ^+X2[7 c, A+, E^])^ , 


where 7 c is the binding momentum in the charged channel, 7 c = \/2 pdd*BEc, where BE^ = 


11 










mr)*± + m£)± — mx{ 3872 ), and go and 5 ^+ are the conplings of the X(3872) to the nentral and 
charged channels. These obey the constraint [H] 

go R'eSQ(—Ex) + 5'+ ReS'^(—Ex) = 1. (12) 


where So(—Ex) and S+(—Ex) are the contribution to the self-energy of the X(3872) from 
the neutral and charged mesons, respectively, and ' denotes differentiation with respect to 


the energy. Eq. (12) can derived by solving the coupled channel problem, see for example 


Ref. |15] where the coupled channel problem is solved for a theory of non-relativistic heavy 
mesons with contact interactions that mediate R-wave scattering in both the 1 = 0 and 
1 = 1 channels. The coupling can be extracted from the residues of the T-matrix at the 


X(3872) pole, which can be shown to satisfy Eq. (12).® If only J = 0 scattering is present 
then go = g+. 

Noting that ReSo(-Ex) = 
can be solved by setting 


— ReS'j_(—Ex) = the constraint in Eq. (| 12 |) 


go = 


2 vr 7 „ 

— cos 0 , gj^ = 
hDD. 


27r7e 


sin^, 


(13) 


SO the decay rates in terms of 9 and the binding momenta are 


r[X(3872) ^ = 


(14) 


‘2g^9t 


2 ^Xco 3 


9irV^'‘°“'mx(3S72) 

r[X(3872) ^ X.i7r“] = 


p® (cos6'y^Eo[7n, Ao, E^] - sin6'y^Eo[7c, A+, E^])' 


9 ‘^gl ,.2 

mx(3872) 


pI (cos6'y^Ei[7„, Ao, EJ - sin6'y^Ei[7c, A+, E^]f 


r[X(3872) ^ Xc27r®] = 


^9^gl 


m 


)Pdd* 


Xc2 ^3 


pi (cos6'y^E2[7n, Ao, E„] - sm 6 y/YcF 2 [lc, A+, E^]y . 


18x2/2— rnxi3872) 

The actual value of 6 depends on the underlying dynamics, and cannot be determined from 
the EFT a priori, so we will leave it as a free parameter. By tuning 9 we can arrange a 
cancellation between charged and neutral loops which allows the prediction to be consistent 
with the bounds. Demanding ^^r[X(3872) —)■ Xcjtt®] < 79 keV, we hnd that 9 = 0.37 ± 
0.04. For this range of 9, 0.78 < 90 / 9 + < 0.99, so the ratio of these couplings is close to 


® I thank R. P. Springer and J. Z. Liu for discussions about this point. 
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FIG. 2: ro/r 2 (solid) and ri/r 2 (dashed) as a function of Q. 


1. This range is consistent with scattering being dominated by the 1 = 0 channel. The 
constraint on 0, and hence 5'o/5'+) is correct so long as 7 „ ~ 14 MeV, and g\ ~ 1/(500 MeV). 
Unfortnnately, the uncertainties on both these parameters are 0(1). If these parameters are 
an order of magnitude smaller, which seems unlikely but is not ruled out by experiment, 
then the constraints on Q and go/g+ would be considerably weaker. 

Finally, we comment on the predicted ratios for Tq : Ti : r 2 in this approach. Because 
the desired rates are achieved by a hne-tuned cancellation between charged and neutral pion 
loops, the ratios vary wildly as a function of 9 near 9 = 0.37 where all three decay rates 
are very close to zero. The plots in Fig. [^show the ratios F 0 /F 2 (solid) and F 1 /F 2 (dashed) 
as a function of 9. The plot on the left in Fig. shows these ratios for a wide range of 9 
and one sees that F 0 /F 2 ~ 3.2 and Fi/Fq ~ 1.2 for most values of 6*, except near 9 = 0.37. 
The plot on the right shows the prediction for the allowed range 0.33 < 9 < 0.41. In this 
range the ratios deviate signihcantly from Fq : Fi : F2::3.2:1.2:l. It would be interesting 
to obtain experimental information on Fq : Fi : F 2 as this could distinguish between the 
various approaches to calculating the X(3872) decays to conventional charmonia. In the 
hadronic loop approach, with both charged and neutral mesons included as explicit degrees 
of freedom, measurement of these ratios could determine the correct value of 9. 


III. FACTORIZATION 

In this section we discuss how the hadronic loop approach discussed in the previous 
section is related to the factorization approach of Ref. [18]. We begin by considering the 
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amplitude from the loop diagram in Fig. [^), ^ with only neutral D mesons in the loop. 
This is given by 


^ib[X(3872) -)■ Xco7r°] = iex ' ^ 


(15) 


dH 


(27r)4 Ex - Ao + /o 


P 

2m*r 


+ ie —/n 


2mD + Ex + /o - ^TT - 


Here Ex is the energy of the X(3872) relative to 2m d, so Ex = Aq The/o 

integral is done by contour integration, resulting in the integral: 

A rW/'QS70^ ^ 0l [^99l f 2/iDD* 1 

Aib[X{3872) XcoTT ] = -ex ■ Pn\ -F- /2, ^2 • (16) 


2 U J (27r)3 P + 7n Ex - ' 


2m D 


Note that the integrand scales as J d^l/{Py for large I and hence the integral is hnite. 
When we multiply this amplitude by the factor y/27r'jn/ Pdd* coming from the wavefunction 
renormalization, this result can be written as 


A1i6[A(3872) —>■ Xco'^^] — ^ ^ ^ib[A(3872) —)■ Xcot^^] 


flDD* 

r dH 


(17) 


(2vr) 


{l)M[D*\l )d\-1)^XcoA, 


where 


) — 


\/8vr7„ 


P + lV 

is the momentum space wavefunction of the D^-D*^ in the A(3872), and 

^99i 1 


(18) 


M[D*\l )D\-l ) ^ = -ex -P. 


2 U Ex-E^- 


d-PnP 


(19) 


2m£) 2m£) 

is a tree-level contribution to the HHyPT amplitude for D^D*^ —>■ XcA- The momentum 
space wavefunction, Pdd*{1 ), has the form dictated by ERT and is correct so long as the D 
mesons are separated by large distances compared to the strong force that binds them. In 
a theory in which the charged D mesons have been integrated out, the scale Xc = 126 MeV 
should be considered large and the wavefunction Pdd*{1 ) can only be trusted below this 
momentum. 


^ The discussion that follows applies to all three diagrams. 


14 


















With this in mind, we will continue evaluating Eq. (16) in the hadronic loops formalism, 
but now imposing a UV cutoff on the integral. Combining the two terms with Feynman 
parameters, we get 


-^i 6[-^(3872) —>■ ] — ex ' P-k\1 2~^ 


^99i 


dx 


= ex -Pr^ 


3 99i 
2 47r2/7 


■2flDD*T^D / dx 


^ d^l 2fi£)£)*Tn£) 
(27r)3 (/2 + A(a;))2 

A 1 


( 20 ) 


+ 


A2 + A(a;) 


: tan 



where A{x) is given by 


Pt 


A{x) = + X moiE^ - Aq) + ^ + 7 . 


2 ~ 2fi£)£)* 


. 2 PI 


X 


( 21 ) 


2 2fiDD* 

The term 'y'^irnD — 2 ixdd*)/{‘2pdd*) ~ —2.8MeV^ is negligible compared to the remaining 
terms so we will drop it as well as the terms proportional to in A(x) since -C mDiE-,, — 
Aq). One can check that setting = 0 only changes the numerical values of the functions in 
Eq. ([^ by a few percent. This approximation allows the integral in Eq. (20) to be evaluated 
analytically and one obtains 


A1i 6[A(3872) —>■ Xco7r°] — ex ■ P 


3 99 i ‘^pDD* 


X 


2 2tt^U E^ - Ao 

y/mD{E^ - Ao) +72 tan”^ 


( 22 ) 


A 


y^moiE^ - Ao) +72^ 

Since the integral is hnite we can send A —)■ cx) and the result is 


7 „ tan I — 
\7n 


Alife[A(3872) —)■ Xco^r ] — ex 'Ptt\ 7 


3 99i 2/i£)D* 


-7n + \/mD{,E^ - Ao) +il . (23) 


2 47r/^ E^- Ao l 

This is the result from the hadronic loops formalism when we set = 0. To see this it 
is helpful to note E{a,b,0) = 2{—y/a + \/a + h)/b. As stated earlier, setting p^ = 0 is 
an excellent approximation to the exact result. However, we have argued that in a theory 
without explicit charged mesons the cutoff A should be not much larger than 7 c ~ 126 MeV. 
The factor mDiE^, — Aq) ~ (736 MeV)^, so ttid^Et, — Aq) 3> 7 c- In Figs. [^) and c), the 
factor of rriD^ET, — Aq) is replaced with m£,(*)(i7^ + Aq) or 2pioD*ET,, which are even larger. 
Typically these quantities are of order (850 MeV)^ and can be as large as (1073 MeV)^. So, 
for a physical value of the cutoff we should take 7 ^ -C A — A), then we have 


Alife[A(3872) —>■ Xco7r°] — ex ■ Pt, 


3 99i ‘^Pdd* 

2 dvr/^ E^ - Ao 


2 A 


TT 


— -7n +0 2^, 


A 


In 


A itidEt, 


.(24) 
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This is actually the result in the factorization approach. Starting with Eq. (16) we note 
that in the second propagator Ex—E.^ — (/^+(/—p,r)^)/(2m£)) Aq—+O( ^ , where 
^ ~ 7n ~ Q (Q denotes any generic scale of order the binding momentum). Consistency of 
XEFT power counting requires that we drop the terms. Then the I integral 


is straightforward and one obtains Eq. (24) for the amplitude. Note that after dropping 


the 0( 


Q ,PttQ,pI 

nioE-n 


) terms the integrand scales as j d?l/l‘^ for large /, hence the integral i 


IS 


divergent and depends on the cutoff. In the factorization formalism, the divergent integral 
is interpreted as the nonperturbative matrix element 

2 


1 


5^|(0|^£j(A)VP + VP)|A'(3872,A)>p = . (25) 


Here the evaluation of this matrix element is sensitive to the cutoff A. This indicates 
the matrix element is sensitive to the short-distance nature of the A(3872) and cannot be 


calculated with XEFT. Still, we can use the formula in Eq. (25) to parametrize the matrix 
element and the constraint on this matrix element from the requirement ^jr[X(3872) —)■ 
XcjE^] < 79 keV, when expressed in terms of A, is A < 325 MeV. This conhrms that in the 
theory with charged mesons integrated out, the cutoff must be interpreted as being a few 
hundred MeV at most and much lower than the scale set by 

In the context of Non-Relativistic QCD, making similar expansions in non-relativistic 
propagators inside loop diagrams in order to maintain consistent power counting is known 
as the multipole expansion |16]. In the present case, this keeps contributions from the loop 
integral that come from low scales I ~ Xn but discards contributions that come from high 
momentum region of integration I ~ A/moE)) ~ 850 MeV. When the cutoff is taken to 
inhnity, contributions from both regions contribute to the hnite answer (see the two terms 
in Eq. (23) ) and the contributions from large I ~ give the dominant contribution. 

It is conceivable that in a theory with explicit charged and neutral D mesons the true cutoff 
can be taken 0{GeV) and this second contribution can be reliably computed. But in a 
theory with only neutral D mesons the cutoff cannot be interpreted as being much higher 
7 c ~ 126 MeV, otherwise the charged mesons should appear as explicit degrees of freedom. 

Finally we consider what happens when the charged mesons are included in the theory. 
Let us assume that in the theory with explicit charged mesons that we can take A to be 
large and keep the region of the integral from I ~ Neglecting terms suppressed by 

^^/(m^^E^), the contribution to the matrix element for A(3872) —)■ Xco^^ from the diagram 
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in Fig. [^), with both neutral and charged mesons and the relevant couplings included, is 


-^ i 6 [-^( 3872 ) —)■ Xco^r^] — ex ' P-n 


3 ggi 

2 47r/^ 


2pDD* 


(26) 


X 


- 7 n + \/mD{E^ - Ao) +7^ 7^ - 

9o --X-+ 9+ - 


Ett — Ar 


Ett — A_ 


In the / = 0 limit, go = g+, the terms proportional to go\/mr,{ET, — Aq) + 7 ^ and 
g+^mD{ET, — A+) + 7 ^ essentially cancel because they differ by only 2% in magnitude. 
Noting that Ao/A+ = 1.01, the final result is well approximated by 


Mib[X{3872) -)■ Xco7r°] — goex ■ Pny 2 4^} E^-^Aq ^ ' 


(27) 


which is the factorization result in a theory with only neutral mesons, Eq. (24), with the 
UV cutoff. A, replaced with A = 7 r 7 c /2 198 MeV. So in this limit the hadronic loops 

result is equal to the factorization result in a theory with only neutral D mesons, with an 
appropriately low value for the UV cutoff. Note that in isospin conserving decays the high 
energy contributions from charged and neutral loops will add rather than cancel. Their 
effects must be reproduced by diagrams with local counterterms in XEFT. 


IV. CONCLUSIONS 

In this paper we studied the decays F [A(3872) —)■ Xcjtt^] within the two commonly used 
approaches to calculating A(3872) decays to conventional quarkonium within EFT: the 
hadronic loop approach and the factorization approach. Within the hadronic loop approach, 
we find that if one only includes neutral mesons as explicit degrees of freedom, and uses 
the estimate of the Xcj coupling to D mesons from Ref. |36], then predictions for each of 
these partial widths exceeds the known bound on the total width. We then obtained a 
bound on ^jF[A(3872) —)■ by exploiting the fact that F[X(3872) —= 

F[£)*o D^tt^] in the limit of small binding energy within ERT. Combining this with 
known results for F [A(3872) —Zi)°i)°7r°]/F[A(3872)] as well as lower bounds on branching 
fractions to observed decays of the A(3872), we found ^jF[A(3872) —)■ Xcj 7 r°] < 79 keV. 
To calculate the theoretical uncertainties in the estimation of F[A(3872) — D^D^n^] we 
use the results of Ref. [9] , which supplement ERT with range corrections, corrections from 
higher dimension operators in XEFT, and pion exchange. We conclude that the prediction 
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for j r[X(3872) —)■ Xcjt^^] is almost two orders of magnitude too large if Xn ~ 14 MeV and 
gi ~ l/(457MeV). Within the hadronic loop approach, the prediction for ^jr[X(3872) —)■ 
XcjT^^] can be made consistent with experiment by including charged charmed mesons in 
addition to neutral charmed mesons as explicit degrees of freedom. The couplings of the 
X(3872) to the neutral (gfo) and charged channels must be tuned to arrange a near 
cancellation between the charged and neutral meson loop contributions. Consistency with 
data requires 0.78 < 5'o/fi'+ < 0.99. If X(3872) appeared as a pole in the J = 0 channel only 
then we would expect (y'o/fi'+ = 1; so the cancellation is naturally explained if the X(3872) 
is an / = 0 state. 

Next we discussed the relationship between the hadronic loop approach and the factoriza¬ 
tion approach to X(3872) decays. We showed that the hadronic loop diagram is proportional 
to the integral /—)■ Xcj7r°], where '^dd*(0 is the momentum 
space wave function of the X(3872) predicted by ERT and —)■ Xcjt^^] is the 

tree-level amplitude for D*^D^ —)■ Xcjtt® in HHyPT. The integrals are well-approximated 
(within ~ 5%) dropping terms that are p\/{rriDET,) suppressed. Making this approxima¬ 
tion, we see that the hadronic loop integral contains two widely separated energy scales: 
Xn = 14 MeV and y/moR^ ~ 850 MeV. The hadronic loop result is numerically dominated 
by large loop momenta of order y/mRE^. For these high momenta, 'ipDD*{l) is likely to 
deviate from ERT form, since this form is only known to be correct for / ~ 7 „. If charged 
mesons have been integrated out of the theory, the ERT form of the wave function is only 
reliable for I < Xc- If this is the case the theory must be interpreted as having a UV cutoff 
~ 100 — 200 MeV. In the limit A the hadronic loop is well approximated 

by the factorization formulae for the decay rate. 

The factorization formulae for the X(3872) decay rate can be interpreted as performing 
the multipole expansion on the hadronic loop integral, i.e., the XEFT power counting / ~ 
In ^ Q ^ V^dE-w is imposed at the level of the integrand in XEFT. Since I ^ Xn ^ Q 
there is no further approximation for 'ipDD*{l), but within Ai[D*^{l)D^{—l) -X Xcj7r°] we 
drop corrections suppressed by {l,Xn)/'m-D and Pn/mo- This has the effect of removing the 
I ~ y/rnoE^ contribution to the hadronic loop and keeping only the low energy I ~ Xn 
contributions. In the hadronic loops approach with both charged and neutral charmed 
meson, if we choose go = g+ so the large / ~ y/mRE^ contributions cancel, the remaining 
terms in the integral are the same as the factorization result in a theory with explicit neutral 
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mesons only, and A = ~ 200 MeV. 

We conclude that within the hadronic loop approach it is inconsistent to keep only neu¬ 
tral charmed mesons and integrate loop momenta to arbitrarily large momentum. If loop 
integrations are taken to inhnity, keeping large contributions from I ~ ~ 850 MeV 

then charged D mesons must be included and the coupling of the X(3872) to the charged 
channel must be nearly equal to that of the neutral channel so the predicted rates for 
^jr[X(3872) —)■ Xcj7r°] are consistent with data. If the charged mesons are integrated out 
of the theory, then the cutoff on the loop momenta should be A ~ 0(%) and the results 
of the hadronic loop approach will be consistent with what is obtained in the factorization 
approach to X(3872) decays to conventional charmonia. 
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